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The purpose of this paper is to analyze in detail the Hamiltonian formulation for 

the compact Gowdy models coupled to massless scalar fields as a necessary first step 

q-i towards their quantization. We will pay special attention to the coupling of matter 

^ | and those features that arise for the S 1 x § 2 and 8 3 topologies that are not present 

in the well studied T 3 case -in particular the polar constraints that come from the 

regularity conditions on the metric. As a byproduct of our analysis we will get an 

alternative understanding, within the Hamiltonian framework, of the appearance of 

CO \ initial and final singularities for these models. 
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I. INTRODUCTION 



■ i 

^ ■ Symmetry reductions are a way to gain useful insights for difficult problems in classical 

and quantum general relativity. In this respect the two Killing vector reductions provided by 
the so called Gowdy models are specially attractive because they have a cosmological inter- 
pretation and share some interesting features with similar reductions such as the Einstein- 
Rosen waves -in particular their solvability- both at the classical and quantum regimes. 
Most of the work on these models, after the initial papers by Gowdy has profusely 

analyzed those corresponding to the 3-torus spatial topology and, in fact, this is by far the 

preferred choice to discuss quantization issues I, i, S u 0,1, S m H 0, H H • The 



other possible closed (compact and without boundary) topologies, the three-handle S 1 x § 2 , 
the three-sphere § 3 , and the lens spaces L(p, q), are interesting in their own right. From the 
physical point of view their most salient feature is the fact that they describe cosmological 
models with both initial and final singularities. For this reason they will become useful test 
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beds for issues related to quantization in cyclic universes. 

The first step towards quantization is the Hamiltonian formulation of the model at hand. 
This is specially so for constrained systems where the identification of the relevant constraints 
is a necessary first step either to attempt a phase space reduction, a gauge fixing, or a 
quantization a la Dirac where the physical Hilbert space is identified as the kernel of suitable 
self-adjoint operators representing the constraints. 

Our goal in this paper is to perform a detailed Hamiltonian analysis of the compact 
Gowdy models coupled to massless scalar fields, extending in several ways previous results 
on this subject [If], 17]. Adding matter fields to the system is a way to enrich these models 
and get closer to physically realistic situations. Here we will work in the spirit of previous 



treatments for other two Killing vector reductions [18|, [19| , paying close attention to the con- 
straint analysis, gauge fixing, and deparameterization. To our knowledge the Hamiltonian 
analysis, for the vacuum S 1 x § 2 and S> 3 Gowdy models, has only been addressed in a partial 
way in [20| where the authors give Hamiltonians for these systems. However they do not 
provide the detailed phase space description (constraints, gauge fixing, and so on) necessary 
to understand relevant geometrical issues. Also their reduced phase space treatment does 
not allow to follow other roads to quantization such as Dirac's approach or the viewpoint 



pioneered by Varadarajan in [21[. Among several issues we want to find out how the topol- 
ogy of the spatial slices affects the definition of the constraints, and how the coupling of 
massless scalar fields is realized in the different topologies. Along the way we also want to 
understand in a detailed way, and within the Hamiltonian setting, the mechanisms leading 
to the appearance of final singularities. Our results can be immediately particularized to 
the vacuum situation. 

The paper is structured as follows. After this introduction we will review in section [TT] 
the main points concerning the Geroch reduction for polarized Gowdy models coupled to 
massless scalar matter fields. In particular we will show that this reduction, and a subsequent 
conformal transformation, allows us to interpret these models as 2+1 gravity coupled to a 
set of massless scalar fields with axial symmetry. Some details related to this reduction vary 
depending on the topologies (in particular those related to the quotient spatial manifold) 
and will be commented separately for each case. 

Section [TTT] will be devoted to discuss the Gowdy T 3 models coupled to massless scalar 
fields extending the previous treatments for the vacuum case 0, [lj] (and similar models 



such as Einstein- Rosen waves 16|, ll3]). Along the way we will clarify some issues related 



to deparameterization and the appearance of singularities in the 3 + 1 dimensional metrics. 
This section will be the basis of the treatment that we will follow to study the other possible 
topologies. 

Section [TV] will be devoted to the Hamiltonian formulation of Gowdy models in S 1 x S 2 
coupled to massless scalars. Here we will have to pay special attention to the identification 
of the regularity conditions that the basic fields describing the model must satisfy as a 
consequence of the regularity conditions on the metric. As we will see the constraints that 
are relevant here are different from the ones present for the 3-torus due to the presence of 
a symmetry axis in the spatial manifold. In particular we will get what we will refer to as 
"polar constraints" involving the values of the basic fields at the poles of the two dimensional 
sphere that appears as the quotient space after performing a Geroch reduction. As we will 
show they are first class and play a relevant role to guarantee the differentiability of the other 
constraints. Another issue that will be discussed is how the deparameterization achieved 
by a partial gauge fixing works for this model and how one can arrive at a reduced phase 



3 



space description. We will see that, as it also happens in the T 3 case, the dynamics of the 
system is described by a time dependent Hamiltonian though the time dependence now is 
different and reflects the appearance of initial and final singularities. In fact, as a result of 
our analysis, we will get a geometric understanding of this fact in terms of the geometry of 
the constraint hypersurface in phase space. 

After this we will perform a similar analysis in section |V] for the § 3 topology. Here the 
main difference stems from the fact that we will be forced to perform the Geroch reduction 
needed to describe the model in 2 + 1 dimensions by using a Killing field whose norm vanishes 
on a circle S 1 . This will introduce some modifications in our description and will change 
the analysis of the relevant regularity conditions for the metric. Nevertheless we will find 
out that the final description is quite similar to the one corresponding to the three handle 
discussed above. 

The detailed quantization of the S 1 x § 2 and § 3 Gowdy models will be carried out else- 
where. A fact that will play a relevant role there is the possibility of describing the compact 
Gowdy models in the different topologies as field theories in certain conformally stationary 
curved backgrounds. As this point of view is also useful to understand some of the issues 
discussed in the paper from a different perspective we will show in section |VT] how this can 
be done. 

We end the paper in section IVHl with a discussion of the main results and suggestions for 
future work on this subject. 



II. GENERAL FEATURES OF COMPACT GOWDY MODELS: GEROCH 
REDUCTION AND 2 + 1 DIMENSIONAL FORMULATION 



Let us consider a smooth, effective, and proper action of the biparametric Lie Group 
:= U(l) x U(l) = {(<7i, g 2 ) = (e lXl ,e tX2 ) | £i,x 2 E R(mod27r)} on a compact, connected, 



and oriented 3-manifold ^E. It can be shown 22|, |23[ that this action is unique up to 
automorphisms of and diffeomorphisms of ®E. The spatial manifold ^E is then 
restricted to have the topology of a three-torus T 3 , a three-handle S 1 x § 2 , the three-sphere 
S 3 , or the lens spaces L(p, q) (that can be studied by imposing discrete symmetries on the 
S 3 case). 

Let us take a four manifold ^Ai diffeomorphic to R x ^ 3 ^E and such that (^Ai, ^g a b) 
is a globally hyperbolic spacetime endowed with a Lorentzian metric 1 ^>g a b- Let us further 
require that G^ acts by isometries on the spatial slices of ^Ai. In this paper we will focus 
on the so called linearly polarized case, hence, the isometry group will be generated by a pair 
of mutually orthogonal, commuting, spacelike, and globally defined hyper surf ace-orthogonal 
Killing vector fields (£ a ,cr a ). 

Let us consider now the Einstein-Klein-Gordon equations 

WiU = 87rGW(d0) a (d0) 6) ^g ab ^V a ^V b ^ = (2.1) 
corresponding to (3+l)-dimensional gravity minimally coupled to a zero rest mass scalar field 



Throughout the paper we will use the Penrose abstract index convention with tangent space indices 
belonging to the beginning of the Latin alphabet 2M . Lorentzian spacetime metrics will have signature 
( — h ++) and the conventions for the curvature tensors are those of Wald 25j . 
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symmetric under the diffeomorphisms generated by the Killing fields = C a (p — 0, 

9db — £j^9ab — 0). Here ^'R a b and ^V ffl denote the Ricci tensor and the metric 
connection associated to ^'g a b> respectively. The exterior derivative of the scalar field is 
denoted by (d0) a and Gn is the Newton constant. 

In order to get a simplified, lower dimensional description we will perform a Geroch 
reduction [26[ by taking advantage of the existence of Killing vector fields. The possibility 
of finding the necessary non-vanishing Killing field £ a will depend, as we will see later, on the 
spatial topology that we consider. In some cases the appropriate Killing vectors vanish on 2- 
dimensional submanifolds but, nevertheless, we will be able to use Geroch's procedure even in 
this situation. The idea is to find a suitable reduction on the manifold ^Ai = ^A4/U(l), 
diffeomorphic to R x ( 2 )£, where ^Ai denotes the set of points in ^Ai in which £ a is 
nonvanishing, and reintroduce the removed points (the symmetry axis) as a boundary where 
the fields must satisfy certain regularity conditions. In the present situation hypersurface 
orthogonality will allow us to view as an embedded submanifold, everywhere orthogonal 
to the closed orbits of £ a , and endowed with the induced metric ^g a b '■= ^9ab — Aj 1 ^^, 

where A 5 := ^g ab ^ b := U a > 0- 

In the linearly polarized case the twist of the Killing fields vanishes and the field equations 
can be written as those corresponding to a set of massless scalar fields coupled to (2+1)- 
gravity by performing the conformal transformation g ab := \^g a b- The system (12.11) is then 
equivalent to 

R "b = g X}( d ^)«( d &)* ' 9 ab V a Vb<f>i = 0, C*g a b = 0, 0^ = 0, (2.2) 

i 

where R a b and V a denote, respectively, the Ricci tensor and the Levi-Civita connection 
associated to g a b (all of them three dimensional objects), we have defined 2 4>\ '■= logA^, 
02 := ^/IQttGn^, and we must remember that we have the additional symmetry generated 
by the remaining Killing vector field a a . Notice that (I2.2p are formally symmetric under the 
exchange of the gravitational and matter scalars. However, it is important to realize that for 
some of the topologies that we will discuss, these fields may be subject to different regularity 
conditions in the gravitational and matter sectors that effectively break the symmetry among 
them. The relevant details for each spatial topology will be given in the corresponding 
section. 

In order to obtain the Hamiltonian formulations for the models that we are considering 
here we want to derive the previous equations from an action principle. To this end we 
introduce the (2+l)-dimensional Einstein-Hilbert action corresponding to gravity coupled 



2 2+1 massless scalar fields will be denoted be the subindex i — 1,2. The subindex i = 1 will label the 
gravitational scalar that encodes the local gravitational degrees of freedom in Gowdy models and the 
subindex i = 2 will label the original 3+1 matter scalar. It is completely straightforward to couple any 
number N of massless scalar fields, in practice this can be done by supposing that the index i runs from 
1 to N. 
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to massless scalars 



167tG 3 1 /( t0jtl ) x (2) E 



/ ^el^l 1 / 2 ^. (2.3) 

o7TCt 3 J{t }x( 2 )su{ti}x( 2 )s 

Here i? denotes the Ricci scalar associated to g a t>. K and h a b are, respectively, the trace of 
the second fundamental form K a b (defined by the exterior normal unit vector n a ), and the 
induced 2-metric on the boundary {t } x' 3 'EU {ti} S. Finally G3 denotes the Newton 
constant per unit length in the direction of the ^-symmetric orbits. We have restricted the 
integration region to an interval [to, ii], where t defines a global coordinate on M. The action 
is written with the help of a fiducial (i.e. non dynamical) volume form ^ 3 ^e compatible with 
the canonical volume form ®e defined by the metric g a b- This is given by 3 ®e = y^L^j^e. 
The volume form ®e induces a 2-form ^e a b = ^t a b c n c on each slice {£} x ( 2 ^E that agrees 
with the volume associated to the 2-metric h a b- We have also introduced a fixed volume 
2-form ®e on {£} x ( 2) £ such that ( 2 )e = y/\h\ ^e, and verifies y/\g~\ ^e abc n c = \f\h\ { ^e a b- 
We require that ®e and ^ 2 ^e be time- independent, i.e. C^e = 0, Ct^e = where C t 
denotes Lie derivative along t a := (d/dt) a . We also demand them to be invariant under 
the action of the remaining Killing vector field. In particular, given the (2+l)-dimensional 
splitting of ( 3 )£ it is natural to choose ®e = dt A ^ 2 ^e, with C t ^e = C a ^e = (25| . 
For all the different topologies, using the Stokes theorem, we get 



{3) S(g ab , fc) = -V / <** [ ! 2)e \9\ 1/2 { {2)R + K abK ab - K 2 - \ V g ab (d<p t ) a (d&) b ) (2-4) 
167rG 3 J to y(2) E ^ 2 J 

where we have used the relation R = + K a bK ab — K 2 + 2V a (n a K — n b Vbn a ) and 
denotes the Ricci scalar associated to h a b- Our strategy in the different topologies 
that we will study in the paper will be to write down an action of this type, adapted to the 
peculiarities of the different spatial topologies (in particular those originating in the different 
sets of regularity conditions that we will have to consider) and use it to derive a Hamiltonian 
formulation for the system. 



III. T 3 GOWDY MODELS COUPLED TO MASSLESS SCALARS 



The T 3 Gowdy model is, by far, the most studied to date both at the classical and 
quantum levels 0, 0, |, |, |, @, 0, S, fi M, Q H H H E We will consider in 
this section the coupling of some types of matter fields and the most important aspects of 
its Hamiltonian treatment paying special attention to the deparameterization and reduced 
Hamiltonian description. 

Let us start by considering the orientable 3-manifold = T 3 = S 1 x S 1 x S 1 , whose 



3 In any basis where the nonvanishing components of ^ 3 ^e have the values ±1, Ig] 1 ^ 2 is equivalent to the 
square of the absolute value of the determinant of the matrix of the metric in that basis. 
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points we parameterize in the form (21,22,23) — (e ie , e ta , e 1 ^) with 9,cr,^ G R(mod27r). In 
particular, we endow ^ 3 ^S with the standard volume form d9 A da A d£. We define the 
following (left) G^-group action 

(91,92) ■ (z!,z 2 ,z 3 ) = (e ix \e^) • (e* e",e*) := (e* e , S xl+ °\ e 1 ^) . 

We can consider now the group orbits defined by the commuting subgroups (91,92) = (e lx , 1), 
( 9l ,g 2 ) = (l,e ix ), xeR(mod27r) 

(e ix , 1) • (e w , e la , e l «) = (e ie , e l{x+a \ e^) , 
(1, e") • (e w , e la , e li ) = (e w , e iCT , e i(iE+0 ) , 

and their corresponding tangent vectors at each point of T 3 obtained by differentiating the 
previous expressions with respect to x at x = 

(0,^,0), 
(0,0,ie*). 

Let us consider the four manifold ^Ai ~ R x T 3 . We introduce now three smooth vector 
fields # a , <r a and £ a , tangent to the embedded submanifolds {t} x T 3 (here t is a global 
coordinate on R). In order to do this, let us fix t G R and define on {t } x T 3 the 3- 
dimensional vector fields 4 (d/d6) a , (d/da) a , and (d/d^) a given in the description of ^Y* at 
the beginning of this section. We extend them to ^Ai by Lie dragging along a smooth 
vector field t a defined 5 as the tangent vector to a smooth congruence of curves transverse to 
the slices {t} x T 3 . Notice that given a one parameter family of diffeomorphisms f t we have 
ft* [£, °~] a — [/t*£, ft*o~\ a so we guarantee that the extended fields commute everywhere. The 4- 
tuple (t a , 9 a ,o~ a , £ a ) defines then a paralelization of ^Ai. Here 9 a is the vector field obtained 
by extending (d/d9) a to the four-dimensional manifold ^Ai; a a and £ a are obtained by the 
same procedure. Once we have introduced these vector fields on ^Ai as background objects 
we restrict ourselves to working with metrics ^g a b satisfying the following conditions: 

1. The action of the group on ^Ai defined by (g±, g 2 ) ■ (t,p) = (t, (g±, g 2 ) -p), t G R, 
p G T 3 , with (g±,g 2 ) • P defined above, is an action by isometries, i.e. £ a and a a are 
Killing vector fields (Cf ] g ab = 0, g ab = 0). 

2. t is a global time function, i.e. ^ g ab (dt)\ ) is a timelike vector field. From now on we 
will consider the manifold ^Ai to be endowed with a time orientation such that this 
vector field is past-directed. 

3. {t} x T 3 are spacelike hypersurfaces for all t G R. In particular := ^g a b£, a £, b > 0, 
K ■= (4 W°* 6 > 0. 

4. £ a and o a are hypersurface orthogonal (this defines the so called linearly polarized 
case). This condition means that the twist of the two fields vanishes. This will 



4 Notice that even though we use a coordinate notation these are globally defined vector fields on the spatial 
manifolds {t } x T 3 . 

5 In particular, take t a := {d/dt) a . 
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ultimately allow us to simplify the field equations and describe the system as a simple 
theory of scalar fields. 

Two simple but important results that can be proved at this point as a consequence of the 
first are the following: 

i) If £ a and a a are Killing vectors and [£, a\ a = then C a (^g ab — £aW\) = 0; 

ii) Furthermore, if we define the vector X a orthogonal to £ a as X a := a a — ^ a {^, b o' b )/X^ it 
satisfies [£, AT] a = and also Cx(^g a b — £a£b/^£) = 0- This means that, without loss of 
generality, we can work with everywhere orthogonal and commuting Killing vector fields £ a 
and a a . In fact, we impose 

5. (9 a , a a , £ a ) are mutually ^g-orthogonal vector fields. 

After we perform the Geroch reduction with respect to the field £ a as described above 
we end up with a set of equations that can be obtained from a 2+1 dimensional action of 
the type (J2~3D with ( 2 )£ = T 2 = S 1 x S 1 . Since the remaining Killing vector field a a is 
still hypersurface orthogonal, and non- vanishing, the corresponding space of orbits ^Ai := 
^A4./U(l) ~ 1 x §' can be identified as an embedded hypersurface in ^ Ai everywhere 
orthogonal to the (closed) orbits of a a . The induced 2-metric of signature ( — h) on ^Ai 
can be written 

s a b = g a b - T~ 2 a a a b , 

where r 2 := g ab o' a o' b = \^\ a > is the area density of the symmetry G^-group orbits. In 
the following we will use the notation r = +V / r 2 . We have now an induced foliation over 
^Ai defined by the global time function t introduced before. Let n a be the g-unit and 
future-directed (g ab n a (dt)b > 0) vector field normal to this foliation, and let 9 a be the g-unit 
spacelike vector field of closed orbits tangent to the slices of constant t, such that 

Q a = e ll2 9 a 

for some extra field 7. If we choose the congruence of curves with t a tangent to ^Ai. Then, 
the congruence is transverse to the foliation, and we can express 

where A^ > and N e are proportional to the lapse and shift functions. The factor e 7,/2 will 
allow us to obtain a proper gauge algebra and simplify later calculations. We require that N, 
N 9 , and 7 are smooth real- valued fields on ( 3 'JA. As we will see in the following the symmetry 
generated by a a will further constraint them, in particular they will be constant along the 
orbits defined by the remaining Killing vector field. The orthonormal basis (n a ,8 a ,<j a /r) 
is positively oriented with respect to the volume 3-form associated to the 3-metric g a b, 
compatible with dt A d# A da, satisfying ^e a b c n a 9 a a a jr = 1. 
The expression of the metric is 

g ab = ei({N 92 - N 2 )(dt) a (dt) b + 2N e (dt) (a (d9) b) + (d6) a {d6) b ) + r 2 (da) a (da) b . (3.2) 
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The fact that the vectors (t a , 9 a , a a ) commute everywhere will translate into necessary con- 
ditions that the vectors n a and 9 a and the scalars N, N , and 7 must satisfy. These are 

C a N = 0, C a N e = 0, £ ct7 = 0, (3.3) 
(C a n) a = 0, (Cj) a = 0, (3.4) 

and 

\{Cei){Nn a + N e 9 a ) - e^ /2 {C t e^ 2 )6 a + Ne l/2 [9,n] a + {C e N)n a + {C e N e )9 a = . 

This last equation can be projected in the directions defined by the basis vectors to give 

l -NC el + C e N + Ne^ 2 n a n b V a 9 b = , 

l -N e C e i + C e N 6 - X -L a + Ne^ 2 9 a 9 b V a n b = , 
6 a a b V a n b = 0. 

These equations are important because they relate the components of the extrinsic curvature 
of some surfaces with derivatives of N, N 6 , and 7. Notice that the scalars <pi are also constant 
on the orbits of a a (the matter scalar 2 because we have imposed this from the start and the 
gravitational scalar X due to the fact that the two Killings £ a and a a commute: C a \^ = 0). 
Therefore, as we will end up with an essentially two dimensional model with fields depending 
only on coordinates t and 9, we will eventually denote C t with a dot and C g with a prime. 
With this convention, we obtain: 

= T~ l e-\r'i - 2t") , (3.5) 
K ab K ab - K 2 = --^(7 - W - 2N e '){r - N e r') , (3.6) 

9 a \d<PU^i) b = - IN 9 Mi + (N d2 - iV 2 )0f ) , (3.7) 

and then the action can be written as 



167rG 



This will be the starting point for the Hamiltonian formalism. Notice that the previous 
expression is coordinate independent. The Lagrangian is written as an integral over the torus 
T 2 of the 2-form obtained by multiplying the fiducial volume form and a scalar function. All 
the terms in this scalar are defined through the use of geometrical objects, in particular the 
derivatives are Lie derivatives along the fields introduced above. This will prove particularly 
important when dealing with other spatial topologies. In this case, it is natural to choose 
as fiducial 2-form (2) e the one verifying ^e ab 9 a a b = iVe 7 r/|5f| 1/2 = 1, i.e. (2) e = d9 A da. 
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The Hamiltonian can be easily obtained by performing a Legendre transform. It has the 
form 

H = C[N]+C [N e ], 

where 6 

C,{N e } = jf N»C e .= 2p' y + Pll '+PrT' + . 

The lapse and the shift act as Lagrange multipliers and enforce the constraints C — 0, 
= 0. The canonical phase space (r, u) is coordinatized by the canonically conjugate 
pairs (7,p 7 ; t, p r ; 0j, and uj denotes the standard (weakly) symplectic form 

u = J ^7 A 5p-y + 5t A 5p T + Hi A <*Pfc) • (3-8) 

The dynamical variables are restricted to belong to a constraint submanifold T c C T globally 
defined by C = 0, Cg = 0. The constraints can be written in an equivalent way by taking 
"linear combinations" obtained by integrating them against suitable weight functions N g and 
Ng in such a way that the vanishing of the weighted form of the constraints for all of them 
is equivalent to the vanishing of C and C$ at every point of S 1 . The gauge transformations 
generated by the (weighted) constraints are 7 

{ 1 ,C[N g ]} = -N g p T) 
{T,C[N g }} = -N gPl , 

{ Pl ,C[N g ]} = -(N g ry, 

{Pr, C[N g ]} = -{Ngi)> + l -N g E - tf) > 

i 

{p* i ,C[N a ]} = (N9TW, 



6 Here and in the following / F := Fd9. 

is 1 Js 1 

7 In the rest of this section we will choose units such that 8G3 = 1. 
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b,C e [N e g ]} = 2N e g ' + N gl ', 
{T,C e [N e g ]} = Ny 7 

{ Pl ,Ce[N 6 g }} = (N d gPi y, 

{p T ,c e [N s g }} = (Nyy, 

{P^C e [N e 9 ]} = (N e gP(t> J. 



A straightforward calculation shows that the constraints are first class in Dirac terminology, 
or equivalently that T c is a coisotropic submanifold of T. Indeed, the Poisson algebra of the 
constraints is a proper Lie algebra 

{C[N g ],C[M g \} = C e [N g M' g - M g N' g ] , 
{C[N g ], Ce[N e g \\ = C[N g M 6 g ' - M e g N' g ] , 
{C 6 [N% Ce[M e g ]} = C e [N e g M e g ' - M e g N e g '} . 

Notice also that, as a consequence of the introduction of the suitable exponential factor e 7//2 
in (13.11) we have a closed gauge algebra [13] (i.e. with structure constants). 

In order to proceed we would like to isolate the true physical degrees of freedom of the 
model. As is well known there are several possible ways to do this. The usual ones are 
gauge fixing, i.e. the isolation of a single point per gauge orbit by imposing appropriate 
extra conditions on the phase space variables, and phase space reduction -that requires us 
to find a way to effectively quotient the phase space by the equivalence relation loosely 
defined as "belonging to the same orbit". The successful implementation of the reduction 
allows us not only to label gauge orbits but also provides us with important mathematical 
structures (topological, symplectic,...) from the ones present in the initial phase space. Here 
we will see that a partial gauge fixing (deparameterization) can provide us with another 
interesting way to deal with the system because it can be described by a time-dependent, 
quadratic, Hamiltonian 0, 0, EH- As we will show below this is also true for the other 
spatial topologies. If one is interested in quantizing the model one can alternatively use the 
Hamiltonian formulation described above to attempt a Dirac quantization. 

The Hamiltonian vector fields associated to the weighted constraints C[N g ], Cg[N g ] are 
tangential to T c and define the degenerate directions of the presymplectic form ui\r c - The 
deparameterization procedure is based on the choice of one of these Hamiltonian vector 
fields to define an evolution vector field Eh r , generated by a reduced Hamiltonian Hr of 
a generically non-autonomous system. With this aim in mind, we will impose gauge fixing 
conditions in such a way that at least one of the first class constraints C is not fixed. This 
will be used to define dynamics. Any remaining constraints left over by the (partial) gauge 
fixing will generate residual gauge symmetries. 

Let l : Yg — > r c denote the embedding of the gauge fixed surface given by the first class 
constraints and the gauge fixing conditions; the pull-back of the presymplectic form to this 
surface, has a single degenerate direction defined by the Hamiltonian vector field Eh r - 
Select then a suitable phase space variable T such that Eh r {T) = 1. The level surfaces of T 
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are all diffeomorphic to a manifold Tr and transverse to Eh r , defining a foliation of with 
T as global time function. In that case, l*uj = — dT AdH R +u R , E Hr = dT+X HR , where u R is 
a weakly non-degenerate form, and (T R , u R , H R {T)) define a non-autonomous Hamiltonian 
system. Any remaining first class constraints will define a constrain submanifold in T R . 



The conditions that are usually considered for this problem 0, E, EE H 



are 



r' = 0, (3.9) 

p; = o. (3.io) 

They mean that both r and p 7 take the same value irrespective of the point of S 1 but they 
do not specify which one. Notice that conditions of the type r = t or p 7 = — p with to,p 
not only would tell us that r and p 7 are independent of 9 but also assign a fixed value to 
them, thus removing additional degrees of freedom. 8 

This means that when using (I3.9|I3.10I) there is still a dynamical mode in r that may vary 
in the evolution -at the end of the day it will be identified with the time parameter- but is 
constant on every spatial slice in the 3+1 decomposition. The fact that this class of models 
have an initial spacetime singularity suggests that there are interesting interpretive issues 
as far as the equivalence of the different choices of gauge fixing is concerned (how does this 
singularity manifest itself after a full gauge fixing? how does it show up if other gauge fixing 
conditions are used?). It should also be pointed out that although it is possible to think of the 
condition r' = as a one parameter family of gauges r — t, with t G (0, oo), it is dangerous 
to use it in this last form when computing Poisson brackets (it would be something like 
"mixing parametric and implicit equations" ) as can be checked by explicit computations. In 
this case the correct attitude would be to work in the extended ("odd-dimensional") phase 
space, mathematically described as a cosymplectic or contact manifold, incorporatin g a time 



variable and employ the usual techniques for non-autonomous Hamiltonian systems [281 . 129 

A convenient way to discuss gauge fixings is to describe our family of gauge conditions 
by introducing an orthonormal basis of weight functions on S 1 

Y n {6) := ^=e m6 , n G Z , 

v27T 



and consider the family of constraints C[y„], C#[Xi]- By expanding now 



with 



-inO _ „ _ 1 / „—in6„ 



T n = -/= / e T > Pin = -7= / e p 7 
V27T Jgi V27T Js 1 



A useful example to appreciate the difference between taking some derivatives to be zero and fixing the 
values of the functions is to consider the straight line X\ = X2 = in R 3 where all the points have equal 
coordinates in contrast with the point x\ = x% = = 1. 
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the previous gauge fixing conditions become 

Tn = , p ln = 



Vn G Z - {0} . 



(3.11) 



In order to see if this is a good gauge fixing (and, alternatively, find out if some gauge 
freedom is left) we compute 



{r n ,C[Y m }} 



/27T 



finmP-yo j { T n, Cg \Y m ] } ~ 



{p 7n ,C[Y m ]}«0, { P , n ,C e [Y m }} 



i n 



OnmP'yQ 



where n G Z — {0}, m G Z, and the symbol ~ denotes equality on the hypersurface defined 
by the gauge fixing conditions and the constraints, the so-called gauge fixing surface Tq C 
T c . Notice that with this way of writing the constraints (without the extra terms that 
would be present if we had not introduced the exponential prefactor in (13. ip ) the gauge 
transformations of r and p 7 only involve these objects themselves. It is convenient to write 
the previous expressions in table form 





T\ — p 71 = 


r_i = p 7 _, = 




C[Y ) 
C e [Y ] 












C[Yi] 


Tk^o 








c[y_x] 

C e [Y_i] 






h^* 













As we can see the only constraints that are not gauge-fixed by the conditions introduced 
above, as long as p 7o ^ 0, are C[l] = and Cg[l] = 0. From now on we will consider the 
sector p 7o < 0. As we can see we have two first class constraints left over by our partial 
gauge fixing 



J si ((2r" - T'r') - w + \ ^ ( P f + rtff) - 

/ ( - % + Pji + Ptt' + Y]^< 

Js 1 v 



(3.12) 
(3.13) 



We can pullback the relevant geometric objects to the submanifold Tq defined by the gauge 
fixing conditions to eliminate some of the variables in our model. Denoting by i : Yq — > T 
the immersion map, the pullback of the (weakly) symplectic form (I3.8P becomes 



d7o A dp 70 + dr A dp T0 + V] / 



Hi A 



(3.14) 
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The pullback of the constraints (13.1 2113. 13p is 

C:=-p 10 Pr + WI (V^ + ^fUo, (3.15) 

C e-=J2 f P*#«°- ( 3 - 16 ) 

Let us look now at the gauge transformations of r generated by (13.151) 

{r ,C} « -p 70 , 
{p 7o ,C}^0. 

If s G (0, oo) parameterizes the gauge orbits we see that on them we have r = ps and 
p l0 = —p, with p > 0. This suggests that a simplification of our model will occur if we 
introduce a canonical transformation where ro and substituted for new canonical 

variables. Indeed, the canonical transformation 



t =TP, p T0 = — , 

7o = ~4=(Q + P fT) , p, = -V^P 



with (Q,P > 0), and (T,px) canonically conjugate pairs, allows us to write 



Finally the canonical transformation (here (Q,P) and (y^p^.) are new canonical pairs [6() 



turns the constraints 




1 

' i 1 



P 



(3.18) 



(3.19) 
(3.20) 



c-p^Ij:/^^)-^ 

and the 2-form (13. 14|) becomes 

l*u = dQ A dP + V / % A <fp w + dT A dp T . (3.21 

, is 1 
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The fact that (13.191) is linear in p T allows us to interpret the 4-tuple ((0, oo) x T R , dt, u> R , H r ) 
as a non-autonomous Hamiltonian system with T = t as the time parameter, restricted to 
verify the global constraint (13.20)) . The reduced phase space T R is coordinatized now by the 
canonical pairs (Q, P; (/p^p^J and is endowed with the (weakly) symplectic form 

u R := dQ A dP + V / dtp, A 5p Vi . (3.22) 

The reduced time-dependent Hamiltonian H R (t) : T R — > R is given by 

1 /" (V% % , + , , 2 



^w-L / f+w). ( 3 - 23 ) 



and the evolution vector field is given by 

dt t 5<fi Yl Sp 



Pi 



This defines the only degenerate direction of (13. 2 1 p . Although the form of the Hamiltonian 
that we have just obtained seems to suggest that the gravitational and matter scalars are not 
coupled, in fact the_ constraint (13.201) shows that this is not the case 9 . Notice also that the 
canonical pair (Q, P) describes a global degree of freedom even though they are constants 
of motion under the dynamics generated by (13.231) . The singularities that must be present 
in this case as a consequence of the Hawking-Penrose theorems j25| can be understood as 
coming from the singular behavior at t = of the Hamiltonian (13.231) . 

Finally, it is possible to recover the original 4-dimensional spacetime from this 3- 
dimensional formulation. First notice that the gauge fixing conditions defining the depa- 
rameterization are preserved under the dynamics if and only if the lapse and shift functions 
iV and N e are constant. By redefining the coordinate 6 as in [9] we can eliminate the shift 
function from the metric. We can proceed in an analogous way for the lapse function to 
make it equal to 1. Once we integrate the Hamiltonian equations corresponding to (13.231) . 
undo the canonical transformation defined above, and solve the constraint Cg = in order 
to obtain the 7 function, we uniquely determine the 3-metric (1 3 . 2 j) . and hence the original 
4- metric. 



IV. S 1 x S 2 GOWDY MODELS COUPLED TO MASSLESS SCALARS 

Let us consider now the three-handle ®S = S 1 x S 2 , parameterized as (e^, e la sm9, cos#) 
with 9 G [0, 7r], £,er G !R(mod27r). Using the group parametrization introduced above we 
can write the G^-group action in the form 

{91,92) ■ (e^e^sin^costf) = {e x \e X2 ) ■ (e l? , e ia sm6, cos9) = {e t(xi+0 , e i{x2+a) sin 9 , cos 6) . 



9 The matter fields act as sources for the gravitational field, hence, the solutions to the Einstein equations 
should depend on the matter content. 



15 



The action of the two U(l) subgroup factors of is 

(1, e ix ) ■ {e^, e ia sin 9, cos 9) = {e^, e i{x+a) sin 9, cos 9) , 
{e ix ,l) ■ {e^,e ia sm9, cos 9) = {e i(x+0 , e ia sin 9, cos 9) . 

The corresponding tangent vectors at each point of ( 3 )£, obtained by differentiating the 
previous expressions with respect to x at x = 0, are 

(0,ie ia sin 9,0), 
(ie* 0,0). 

As we can see the second one is never zero but the first one vanishes at the poles of the 
sphere § 2 where 9 = 0, 7T. This corresponds to the circumferences given by (e l «,0, 1) and 
(e^,0, —1). It is straightforward to verify that both fields commute. In view of all this we 
perform a Geroch reduction by using the non-vanishing Killing. After a suitable conformal 
transformation the field equations can be derived from an action of the form (12.41) with 
( 2 )£ = § 2 . All the fields in this action are defined on S> 2 and are symmetric under the 
symmetries generated by the remaining Killing a a . Since this Killing vector vanishes at the 
poles of the sphere S> 2 we cannot build an everywhere orthonormal basis that involves this 
vector. In fact, we know that as S> 2 is not paralelizable this is impossible on general grounds. 
We nevertheless will consider the triplet of vectors (n a , 9 a , a a jr) whenever it is different from 
zero (for all 9 ^ 0, 7r). Taking again the definition fl3.ll) . the form of the metric is the same 
as in the T 3 case (13. 2p . The symmetry of the problem implies also that (N, N e , 7, r, (pi) are 
constant on the orbits of the Killing field a a . 

A very important issue now is the regularity of the metric. From a classical point of 
view the final outcome of the Hamiltonian analysis of the system is a set of equations whose 
solutions allow us to reconstruct a four dimensional spacetime metric and a set of scalar 
fields satisfying the coupled Einstein-Klein Gordon equations. This means that once we 
decide the functional space to which this metric belongs this will imply that the objects 
that appear during the dimensional reduction, gauge fixing and so on may be subject to 
some regularity conditions. In the T 3 case these are simple smoothness requirements but 
in the present case, due to the existence of a symmetry axis, these are more complicated. 
The regularity conditions that the metric components for an axially symmetric metric must 



verify can be deduced as in [23|, |27|. By using the coordinates (t,9,a,^), we can write the 
original 4-metric ^g a b as 

{4) g ab = e^[(N 92 - N 2 )(dt) a (dt) b + 2N 9 (dt) {a (d9) b) + (d9) a (d9) b ] (4.1) 
+ r 2 e-^(da) a (da) b + e*(d£) (d£) 6 ■ 

This means that we have the following regularity conditions for 9 — » 0, tt (if we impose 
analyticity, otherwise we need only to know the asymptotic behavior for small values of 
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sin 9) 

e (7-<M (^2 _ at2) = A ( ti cos Q) j (4.2) 
e (7-0i) iV e = 5 ( t> cos ) sin q ; (4_ 3 ) 

e^ 1 = C(i,cosfl) , (4.4) 
e 7-*i = D(t, cos 0) + E(t, cos 0) sin 2 , (4.5) 
rV 01 = sin 2 6[D(t, cos 0) - £(t, cos 9) sin 2 0] , (4.6) 

where A, B, C, D, E are analytic in their arguments (despite of the fact that they also 
depend on t, as we will use them in the Hamiltonian formulation of the model we will not 
write the t dependence explicitly in the following). Notice, in particular, that the functions 
D and E both appear in the last two expressions. The conditions for the fields themselves 
(dropping the t dependence) are 

<j>i = ii(coa9), (4.7) 

7 = 7 (cos#), (4.8) 

N 9 = N e (cos 6) sinO, (4.9) 

N = N(cos9), (4.10) 

r = f (cos#)sin#, (4.11) 

, , D (cos 6) - E( cos 6) sin 2 8 ~ . . 

r e" 7 = — ) 5- sin 2 , 4. 12 

D (cos 9) + E (cos 9) sin 2 9 ' V ' 

where 0,, 7, iV 61 , N, T : [—1,1] — > K (iV > 0) and can be written as functions of 
A, B, C, D, E. They must be differentiable in (—1, 1) and their right and left derivatives at 
±1 must be defined (equivalently they must be C°° in (—1, 1) with bounded derivative). Sev- 
eral comments are in order now. First we have been able to write all the relevant fields in such 
a way that their singular dependence has been factored out (sin 9 is not a smooth function on 
the sphere) . The functions defined on S 2 as 0, o cos #,70 cos 9, N e o cos 9, N o cos 9, T o cos 9 
can be alternatively viewed as analytic functions on the sphere invariant under rotations 
around its symmetry axis and can be considered as the basic fields to describe our sys- 
tem. In fact we will do so in the following. We will refer to these functions on the sphere 
as (f>i, 7, N , N, T (without the ocos# that will only be used if the possibility of confusion 
arises) and collectively as the hat-fields. In the following we will write everything in terms 
of them. Second we can see that condition (14.121) implies that 

f (±1) = e 7(±1)/2 . (4.13) 

This means that the values of the fields T and 7 at the poles of the sphere are not independent 
of each other. This is a new feature, not present for the T 3 topology, that must be taken 
into account. As we will see these are necessary ingredients to guarantee the consistency of 
the model. 

Given a smooth (and axially symmetric) function on S 2 its Lie derivative along 9 a , Cg, 
cannot necessarily be extended as a smooth function on the sphere. The function cos 9 itself 
is an example of this because Lq cos 9 = — sin 9. We can, however, define a smooth derivative 
/' for a smooth axially symmetric function as the extension of /' := —-r^dof to S 2 (this is 
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formally done by considering / as a function of cos 8 and differentiating). In the following 
the prime symbol will always refer to this derivative. 

It is natural to consider in this case ^ 2 ^e as the fiducial 2- form associated to a round 
metric, such that {2) e ab 6 a a b = Ne*f sin 6/\g\ 1/2 = sin#, i.e. (2) e = sin0d0 A da. Taking this 
into account we get an action 



167rG, 



dt I V>e 



3 Jt JS 2 



iV[(7 f ' - 2T") sin 2 9 + (6T ' - ff ) cos# + 2T] 



[N e f cos - f - A> 9 T ' sin 2 6} [7 + (2JV 0/ + N e *f) sin 2 - 2iV e cos 6} 



2N ^ 



+ 2N°<j>i & sin 2 + (N e2 sin 2 - N 2 )(j>' 2 sin 2 



As we can see it is expressed as the integral of a smooth function on the sphere. This is 
so because all the fields that appear in the integrand are either the hat-fields, their prime 
derivatives or smooth functions of cos 9. The Hamiltonian can be readily derived from the 
previous action and, as in the T 3 case, is of the form H = C[N] + C e [N e ] with 



C[N] = f {2) eNC 
Js 2 



Is 



^eN 



16irG 3 pjPf + 



1 



167TG3 o 

— 5— Pi. 



167rG 3 
f 



(4.14) 

[{2T" - 7 f ') sin 2 6 + (ff - 6T ') cos - 2T] 
^' 2 sin 2 



C e [N 6 



Is 



+ 2^\ / U>vnV, 
{2) eN e C e (4.15) 

( 2 )eiV e ([2p^ - fp-y - ^'Pf - E <%J sin2 + P>t - 2 Pt] co « ^ ■ 



Again, the dynamical variables are restricted to belong to a constraint surface T c C T in 
the canonical phase space of the system (T,ou), globally defined by the constraints C — 0, 
Cg = 0. T is coordinatized by the conjugated pairs (7 SP^T \Pf',(f>i,P^.) and endowed with 
the standard (weakly) symplectic form 



y (2) e ^7 A % + 5T A %, + ^ 50i A 



(4.16) 
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The gauge transformations generated by these constraints are 10 

tf,C[N g ]} = -N gPf , 
{f,C[N g }} = -N gP ^, 

{$ i: C[N g ]} = ^p $i , 

{p^ C[N g }} = N' g {f cos 9 - T sin 2 9) + N g (f + 3T ' cos 9 - f" sin 2 0), 
{p f ,C[N g }} = A>(2cos0-f sin 2 6) + iV g (f cos ^ - 7" sin 2 0) - 2N g sin 2 6 

{P&, C[N g ]} = N' g Tfo sin 2 9 + N g [(f'ti + ffy sin 2 9 - 2T f> cos 9] , 



and 



-2N e g sin 2 9 + N° g {2 cos 9 - f sin 5 
7V e (f cos^-f'sin 2 ^), 



{7,0^]} 

{r,c e [^]} 

{^C e [N e g ]} = -N e g ^sin 2 e, 
{p^ C e [N g ]} = N e g (2p^ cos 9 - p'^ sin 2 9 
{p f , C e [N 9 ] } = N 9 (p f cos 9 - p' f sin 2 9] 
{Pfr C e [N 6 g ]} = N e g {2p k cos9 - p' k sin' 



- N 9 'p^ sin 2 9 , 
N 9 'p f sin 2 9, 



-N 



Bl. 



sin 2 9 



As is the T 3 case, T c C T is a first class submanifold as can be seen by computing the Poisson 
algebra of the (weighted) constraints 

{C[N a ],C[M g ]} = C e [M g N' g - N g M' g ] , 

{C[N g },C e [M 9 }} = C[(M 9 N' g - N g M 9/ ) sin 2 9 + N g M 9 cos 9} , 
{C e [N% C e [M°}} = C e [(M d g N 9 ' - N d g M 9 ') sin 2 9) . 

We want to check now the stability of the "polar constraints" (Te~^ //2 )(±1) = 1. To this 
end we compute 

{fe~^ 2 , C[N g }} = 87rG 3 N g e-^ 2 (fp f - 2p^j , 
{fe-^ 2 , C e [N e g ]} = e-^ 2 [TN 9 ' + A^Tf - f ')] sin 2 9 . 

The first expression vanishes at the poles as a consequence of the constraint (14.151) for 
9 = 0, ri (sin9 = and |cos#| = 1) whereas the second vanishes because of the sin 2 9 



In the following l6irGs = 1. 
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factor. We then conclude that there are no secondary constraints coming from the stability 
of the polar constraints. An interesting point to highlight here is the fact that these polar 
constraints are necessary conditions for the differentiability of the constraints (I4.14f4.15l) . 

Deparameterization in this case is carried out by basically following the same steps as in 
the T 3 case. Again, in view of the gauge transformations, we begin by choosing gauge fixing 
conditions similar to (13.9113. 101) 



f' = 0, (4.17) 
P~=0. (4.18) 

We introduce now an orthonormal basis of functions on the subspace of axially symmetric 
functions on S 2 

2n + l^ 1/2 



Y n (6) = [-^) Pn(cos9) , iigNU {0} 
where P n are the Legendre polynomials. By expanding now 

oo oo 

T = T n Y n , p/y = ^ V^irXn i 

n=0 n=0 




(%P n (cos6)f, P% = (j^^j ' ^ 2 (2) eP n (cos%^, 



f n = = P% , M n E N . (4.19) 



In order to see if this is a good gauge fixing (and, alternatively, find out if some gauge 
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freedom is left) we compute 



{f n ,C[Y m ]}i 
{f n ,C e [Y m }} 



(2n + l)(2m+ 1) 



(47T) 



Is 2 



(2) eP„(cos60-P™(cosf 



Arc 



r.. /(*> + D(*» + 1) im ecosePn(cose)Pm(cme) 



47T 

(m + l)f 



^4^(2™ + l)(2m + 3) 
mPo 

v /4vr(2m + l)(2m- 1) 
otherwise 



if n = m + 1 
if n = m — 1 



{p 7n ,C[Y m ]} a T ^/ (2n + ( 1 4 ) ^ 3 m+1) / ( 2 )eP n (cos£)[P m (cosfl) + costfP/Jcos. 



s 2 

if m = or m < n 
-T - 



'A-n 

-k otherwise 



if m = n 



\r ) r \_' n ; 1 l^eP n (cos9)[P m (cos9) - sin 2 flP^cosf 



(4tt)3 j§2 

ifm = 0orm<n-l 
* otherwise 



where the symbol ~ denotes that we are restricting ourselves to points in the hypersurface 
Tq C r c defined by the gauge fixing conditions and the constraints. The * and * symbols 
denote terms (computable in closed form but with somewhat complicated expressions) that 
are not needed in the following discussion. As before it helps to display the previous result 
in table form 





T\ = p 7>1 = 


T 2 = p 7i2 = 




C[Y ] 
C e [Y ] 




To n 

2^ U 








C[Yi] 
C e [Y x ] 


P-yO T 

* 




-^r 

V Ion 




cm 

C e [Y 2 ] 


* 

t°— * 


Pio 3T 
2^ 

* 













One must also check if the polar constraints are gauge fixed by our conditions (14.191) . To 
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this end we compute 

{f n ,fe-^ 2 }^0, 



The last Poisson bracket is different from zero at the poles (9 — 0, 7r) for all values of n G N. 
As we can see the only constraint that is not gauge-fixed by the conditions introduced above, 
as long as p 7o ^ and T 7^ 0, is C[l\. This is different from the situation in the T 2 case 
where we were left with two constraints instead of just one. 

As we did before we can pullback everything to the phase space hypersurface defined by 
the gauge fixing conditions. The induced 2-form becomes 

l*uj = d% A dp % + df A dp fo +J2 (2) eS(j)i A 5 P<Pi (4.20) 
and the remaining constraint is 

C ■= -P70PT0 + To [4>/5F( log - 1) - 7o] (4.21) 

V 47T 



The gauge transformations generated by this constraint in the variables T and are 

{f ,C} = -p^ , 

0*70 > C } = T 0, 

so if we parameterize the gauge orbits as before with s G (0, ir) we find now T = psin s and 
= — pcoss, p 7^ 0. In the spirit of the previous section we introduce now the following 
canonical transformation ((Q,P) and (T,p T ) denote canonically conjugate pairs) 

T = P sin T , pf o = ^ cos T — Q sin T , 

Pt 

lo — —Q c °s T — — sin T , p^, = — P cos T . 

In addition, as we did in the T 3 case, it is possible to write the remaining constraint in a 
more pleasant form by performing a further canonical transformation (here, again, (Q, P) 
and (ipijPp.) are canonical pairs) 

Q:=PQ + iV / (2) e^0 l; P:=logP, 
Z i Js2 

tpt = tAxr^VPh , p w = {^f /iE j= , (4.22) 
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giving 



C = p T + 4^ e p (log S -^L + P-l)smT + -Y [ <-% ( + sinT sin 2 o] « . 

V47T 2^J S 2 \smT J 

(4.23) 

It is now obvious the interpretation of the system as a non-autonomous Hamiltonian system 
((0, 7r) x T R ,dt,u R ,H R ), where T R denotes the reduced phase space coordinatized by the 
canonical pairs (Q,P',ipi,p Vi ), endowed with the standard (weakly) simplectic form (13.22p . 
The dynamics is given by the time dependent Hamiltonian H R (t) : T R — > M 

H R (t) =4v/77/(log^£ + P-l)sint + ^ Y [ W e (^ + ri 2 smtsm 2 e) , (4.24) 
V47T 2 i Js 2 \smt J 

with the evolution vector field 

^ = - + 4^(log^ + P)sint^ 

+ E / (2)e (-+1T- + (sin 2 ^)'sint-M . (4.25) 

Several comments are in order at this point. First we can see now that the final description of 
our system is somewhat simpler that in the T 3 case because we do not have any remaining 
constraints and the fields <pi and ip 2 are decoupled (at variance with the previous case). 
On the other hand we see now that the dynamics of the global modes, though easy to 
get in explicit form, is not as simple as the one found for the torus. Notice also that the 
Hamiltonian (14.241) is singular whenever sint = 0. This means that if we pick the initial 
time to G (0, ir) in order to write the Cauchy data we meet a past singularity at t = and 
a future singularity at t = n. 



V. S 3 GOWDY MODELS COUPLED TO MASSLESS SCALARS 

Let us finally consider the case where the spatial slices have the topology of a 3- 
sphere § 3 , described as S 3 = {(zi,z 2 ) G C 2 : \zi\ 2 + \z 2 \ 2 = 1}. Let us define the following 
action of on ®S 

(9i,92) ■ (zi,*) = (e lx \e^) ■ (z u z 2 ) = (e^ Zl ,e^z 2 ). (5.1) 

The action of the two U(l) subgroup factors is 

(e ix ,l)-( Zl ,z 2 ) = (e ix z u z 2 ) 
(l,e ix )-(z 1 ,z 2 ) = (z 1 ,e ix z 2 ). 
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The corresponding tangent vectors at each point of S 3 , obtained by differentiating the pre- 
vious expressions with respect to x at x = 0, are now 

(i*i,0) 
(0,iz 2 ). 

As we can see they vanish at z\ = and 22 = (i.e. at the circumferences (0, e^) and (e lCT , 0), 
£, a G M(mod27r)). This fact poses now the question of how one can possibly use them to 
perform a Geroch reduction that requires us to have at least a non- vanishing Killing vector 
field. On some other respects they present no problems, in particular they are commuting 
fields. A useful parametrization of § 3 is z\ = e l(T sin(6 l /2), z 2 = e*^cos(0/2) with 9 G [0, 7r], 
£, cr G M(mod27r), with the commuting Killing fields <r a and £ a given by cr a = (d/da) a and 
£ a = (<9/<9£) a . This allows us to view the three-sphere as a filled torus in which the points 
on the same parallel of the surface are identified (so that the surface itself can be viewed as 
a circle S 1 ). This is helpful to perform the Geroch reduction. 

The fact that the Killing vectors that we have chosen vanish alternatively in two different 
circles poses a problem as far as the Geroch reduction is concerned because to perform it 
one should use a non-vanishing vector. We will show now that the the fact that £ a only 
vanishes in a one dimensional submanifold will effectively allow us to use them to carry out 
this reduction. To this end we start form an action in four dimensions defined on ^Ai, 
topologically Ix S 3 , and remove the circle where the Killing vanishes from the integration 
region. As this is a zero-measure set the integral will not change. Of course one must take 
now into account the fact that the fields in the new integration region cannot be completely 
arbitrary but should be subject to some restrictions (regularity conditions) reflecting the 
fact that they should extend to the full ^Ai in a smooth way. 

Topologically the two dimensional manifold that appears in the action is (2) S = D, 
where D denotes an open disk. The regularity conditions on the disk boundary are such 
that the fields (of any tensor type) behave in a "radial coordinate" 9 exactly as an axially 
symmetric field would do in the axis. Eventually this will allow us to change the disk by a 
two sphere. 

As in the previous cases we are going to use (t a , 9 a , a a ) as coordinate vector fields. We 
will write now 9 a = f9 a , t a = f(Nn a + N e 9 a ) with iV > and / > 0. The scalars /, N, 
and N e are supposed to be smooth fields oeRxD subject to some regularity conditions 
that will be specified later. Notice that we write now / instead of e 7 /2 as in previous cases 
because we want to allow / to go to zero at the disk boundary. The same argument that 
we used for the two previous cases tells us now that 

NC e f + fCoN + Nf 2 n a n b V a 9 b = , (5.2) 
N e C e f + fC e N 9 - C t f + Nf 2 9 a 9 b V a n b = , (5.3) 
9 a a b V a n b = . (5.4) 

The form of the 3-metric g ab is basically the same as in the other cases 

9ab = f 2 [(N 92 - N 2 )(dt) a (dt) b + 2N\dt) {a (d9) b) + (d9) a (d0)b] + r 2 (da) a (da) b 

and the determinant is now given by \g\ = r 2 A^ 2 / 4 . Again (N, N e , 7, f ), are constant on 
the orbits of the remaining Killing field a a and hence they only depend on the coordinates 
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(t,9). Using the coordinates system (t,0,a,£) we write the original 4- metric ^g a b as 

(4) 9 ab = A(iV e2 - N 2 )(dt) a (dt) b + 2N\dt) (a (dO) b ) + (d9) a (d0) b ] (5.5) 
A 5 

A ? 

We have to find out the regularity conditions satisfied by this metric. At 9 = the regularity 
conditions should be of the same type as the ones that we have already used in the S 1 x § 2 
case. Here, however, we also have to impose regularity conditions when we approach the 
boundary of the filled torus that we obtained by removing the circle where the Killing used 
to perform the Geroch reduction vanishes. This can be formally achieved by changing 9 for 
n — 9. By doing this we find 

^(N 02 -N 2 ) = A(t, cos 9), (5.6) 

A ? 

f 2 

?-N d = B(t,cos9)sm9, (5.7) 
A £ 

A f = 4 cos 2 (#/2) [F(t, cos 9) - G(t, cos 9) cos 2 (#/2)] , (5.8) 
f 2 

J — = D(t, cos 9) + E(t, cos 9) sin 2 (#/2) = 

F(t, cos 9) + G(t, cos 9) cos 2 (#/2) , (5.9) 

— = 4sin 2 (#/2)[ J D(t,cos#) - £(t, cosfl) sin 2 (#/2)] . (5.10) 

Here the functions A, B, D, E, F, and G are analytic in their arguments. Notice that they 
are not independent because they are constrained to satisfy 

D(t, cos 9) + E(t, cos 9) sin 2 (#/2) = F(t, cos 9) + G(t, cos 9) cos 2 ''(9/2) . 

We have used the functions sin(#/2) and cos(^/2) because they alternatively vanish on the 
circles where the Killings themselves become zero and have the dependence of a regular 
scalar function in terms of the "radial" coordinates 9 or ix — 9 on the circles where they do 
not vanish. The cosine dependence of the other functions is dictated by regularity at the 
two circles. This is very important because we will be able to write down our model in terms 
of them and, having cos# as their argument they can be interpreted as functions on § 2 as 
in the S 1 x S 2 case. The conditions that the fields must satisfy (dropping the t-dependence) 
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are 

A 5 = e ^ l(cose) cos 2 (6>/2) , (5.11) 

2 = 2 (cos#), (5.12) 

/ = cos(#/2)e^ (cos0)/2 , (5.13) 

N 9 = TV (cos 9) sin 9 , (5.14) 

N = N(cos9), (5.15) 

r = T (cos 0) sin 0, (5.16) 
D(cos9) - E(cos9) sin 2 (#/2) 
L>(cos#) + E(cos9) sin 2 (#/2) 



2|l _- F(cosfl)-G(cosg)cos 2 (g/2) 

F(cos#) + G(cosfl)cos 2 (#/2) ' 1 J 

where we have used sin# = 2 sin(#/2) cos(#/2). Here, as in the S 1 x S 2 case, we have that 
0i,7,iV 9 ,iV,f : [-1, 1] -> R (N > 0). They can be written as functions of A, B, D, E, F, 
and G. They must be different iable in (—1, 1), and their right and left derivatives at ±1 must 
be defined (equivalently they must be C°° in (—1, 1) with bounded derivative). Conditions 
( 15TT71I5TT81) imply that 

f (+i) e -^ +1 )/ 2 = 1 and e ^i(-D-7(-D = 4 . 

These are the polar constraints for the § 3 topology. This is slightly different from previous 
examples because now these conditions involve different pairs of objets at the two poles of 
§ 2 . 

Our starting point is now the action 

(3) %a6,<&) = -V f\t [ Ve\g\^(wR + K ab K ab -K*- V$>^) a (d^ 
167rG 3 Jt \ 2 i 

with <pi = logA^ — 4>x + log cos 2 (9/2) . Notice that we have changed the integration region 
to S 2 because, as we will see, it can be written in terms of the hat-fields that are smoothly 
extended to § 2 . 

As in the case of the three-handle we choose the fiducial volume element to be 
compatible with the auxiliary round metric on the 2-sphere S 2 , i.e. ^ 2 ^e = sin 6^0 A da, with 
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{2) e ab 9 a a b = Npr/lgl 1 / 2 = sin#. In terms of the fields (N, N e , 7, f, 0j) the action becomes 

( 3 ^(iV,iV f? , 7 ,f,0 J ) = 

1— [dt [ (2) ef A>[( 7 'T" - 2T") sin 2 9 + (5T ' - 7 'T ) costf + f ' + -f ] 



167tG 3 

+4[iV e f cosfl-T -iV^T' sin 2 0] [7 + (2N 01 + JV^-y') sin 2 # + (1 — 3 cos 9)N e ] (5.19) 
2iV V 



i 2 + 2N 9 (j) i # sin 2 + (iV 6 * 2 sin 2 9 - N 2 )^ 2 sin 2 



+i[2(l - cos0)(JV^i + (iV e2 sin 2 - iV 2 )^) + (I — cosfl) 2 iV e2 | 

It is important to remark at this point that the action is the integral of a smooth function 
on the sphere. We arrive at this result after several non-trivial cancelations of terms that 
would diverge at the poles. This reflects the fact that indeed, by removing the circle where 
the Killing vector field used in the Geroch reduction vanishes, we arrive at a consistent 
description of the model. It is also worthwhile pointing out that the structure of the action 
is very similar to the one found in the S 1 x § 2 case but not exactly the same, in fact we will 
see later that the differences are important to guarantee, for example, the stability of the 
polar constraints in this case. 

The Hamiltonian of the system can be readily obtained. As in previous cases it can be 
written as a sum of constraints H = C[N] + Ce[N e ] with 



C[N] = / (2) eA>C (5.20) 
is 2 

= f (2) eA r (-WnGsPWf + — V P?" - i'f') sin 2 9 + (ff - 5T ') cos 9 - \f - f] 
J § 2 \ 167rG 3 2 

1 ^ /167TG3J9 2 ; f \ \ 

C e [N e } = [ WeN e C e (5.21) 
is 2 

= jf (2) e#* {[% - 1% ~ f' Pf - foii s^ 2 9 + [fpf - ^ + p A ] cos 9-p^-p^. 

The two previous expressions, together with the conditions at the poles T(+l)e~^ +1 )/ 2 = 1 
and e 2 ^ -1 ) - ^ -1 ) = 4, define the constraints of the system. As before, the polar constraints 
are necessary conditions to guarantee the differentiability of the (weighted) constraints C[N g ] 
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and Cg[Ng]. The gauge transformations defined by C[N g ] and Cg[N g ] are 11 

{*f,C[N g ]} = -N gPf , 
{f,C[N g ]} = -N g p^, 

{kC[N g ]} = N g P -^, 

{p^C[N g }} = N' g (fcos9-f'sm 2 9)+N g (3f'cos9 + f-f"sm 2 9), 

{p f , C[N g ]} = N g [-- 4>[ + (f + fa) cos 9 - 7" sin 2 6] + N' g (3 cos 9 - 1 - 7' sin 2 

-2^sin 2 ^ + ^^(^-sin 2 ^ 2 ) ; 
{P^C[N g }} = [N g f(fa 2 sin 2 + 1 - cos 6)]' , 



fc 2 ,C[iV g ]} = (iV s r0' 2 si n 



2 



and 



{7, C^JVj]} = -2JVJ' sin 2 9 + N e g (3 cos - f sin 2 9-1), 
{f , C e [N d g ] } = iVj (f cos - f" sin 2 9) , 
{0i , Cg [N g ] } = N e g (cos 9-1 -fa sin 2 0) , 
{fa,C [N e g ]} = -N 9 g fa 2 sm 2 9, 
{p^C e [N e g ]} = -(N e gP ^ sin 2 9)', 

{p f , C e [N e g }} = N e g ( Pf cos 9 - p' f sin 2 9) - N e g 'p f sin 2 9 , 
{PL,Ce[N*\} = -{Nt 



The Poisson brackets of these constraints give exactly the same result that we obtained 
for the S 1 x §> 2 topology and, hence, define a fist class constrained surface r c C T. Here 
(r, to) denotes the canonical phase space of the system, coordinatized by the canonical pairs 
(7,^; T,pf, fa, pi.), and endowed with the standard (weakly) symplectic form (14.161) . We 
must check now the stability of the polar constraints. We do this by computing 

{f e^' 2 , C[N g }} = ±N g e-V 2 (fp f - 2pJ , (5.22) 
{f e"^ 2 , C e [N e g }} = e~^ 2 (^N e g f(l - cos 9) + (N 9 g 'f - N e g f + ^Tf) sin 2 d\ , (5.23) 

N9 

{e 2 ^~\ C[N g }} = -je 2 ^{2p k + Tp f ) , (5.24) 
{e 2 ^~\ C e [N 9 g }} = e 2< ^( - N e g (l + cos9) + (2N* - 2N e g fa + N e g Y) sin 2 6) . (5.25) 



11 Again we take 1&itGz = 1. 
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The constraint (15.211) at the poles 9 = 0, it gives respectively, T(+l)pf(+l) — 2p^(+l) = 0, 
and T(— l)pf(— l)+2p^(— 1) = 0. These guarantee that the Poisson bracket (15.221) . vanishes 
at 9 = and (I5.24p vanishes at 9 = ir. The vanishing of (15.231) at 9 = is due to the presence 
of the factors 1 — cos 9 and sin 2 9 and, finally, (15.251) is zero at 9 = it due to the factors 1+cos 9 
and sin 2 9. As in the S 1 x S 2 we conclude that there are no secondary constraints coming 
from the 

stability of these polar constraints. 

The deparameterization in this case follows closely the one for S 1 x S 2 . The same gauge 
fixing conditions work in our case now. The only new element now is checking if the polar 
constraints are gauge fixed or not and this only requires the computation of 



= e 2 ^-y^lp n (cos#) 

which is different from zero at the poles. As we see the situation now is completely analogous 
to the previous case. The pull-back of the symplectic form to the phase space hypersurface 
defined by the gauge fixing conditions is given again by (14.201) . We are left only with the 
constraint 

C := - p%Pfo + f (^\og^=-%- v^(21o g 2 + 3)+0 lo ) (5.26) 
1 



+ £ 




The gauge transformations generated by this constraint on the variables T and p^ Q are the 
same as for the three-handle and, hence, we can use the canonical transformations introduced 
at the end of the previous section to rewrite (15.261) as 

rri q 

p T + (47r) 1/4 e p/2 ^i sinT + 2v / ^e^(log^ + P-log2- J)sinT (5.27) 

V47T 2 

-2 



The description of the system by a time- dependent Hamiltonian is now straightforward. It 
is interesting at this point to compare the dynamics of this model and the S 1 x S 2 one. First 
of all we see that the global mode have a different behavior now, in particular couples to 
cpi through the term e p l 2 Lf>i sinT in (I5.27p . As we see the gravitational and matter modes 
cease to play a symmetric role in this particular description, at variance with the other 
topologies. However, by writing the regularity condition (15. lip with an extra T (as will be 
justified in the next section) it is possible to restore the symmetry between the gravitational 
and matter scalars in a straightforward way. 

As in the previous cases, it is possible to interpret the system as a non-autonomous 
Hamiltonian system ((0,7r) x r#, dt, ujr, Hr), where Tr denotes the reduced phase space 
coordinatized by the canonical pairs (Q, P; fijPipi), endowed with the standard (weakly) 
simplectic form (13.221) . The dynamics is given by the time dependent Hamiltonian Hn(t) : 
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r, 



H R (t) = (47r) 1/4 e p/ VioSint + 2v^Fe p (log^^ + P-log2- -)sint 



sint 



47T 



(2) e f 



+ y?f sin i sin z 9 



sint 



with the evolution vector field given by 



Eh r = — 



d_ 

di 



(47r) 1/4 e p/ Vi sint + 2v^e p sint log -= + P - log 2 - - 



sint 



(47r) l/4 e P/2 sint 



d 



E 



<-" r f + ( S in 2 ^)'sint— 



sint <5(^j 



5p, 



<Pi 



_d_ 

dQ 



The singularities in this case show up in the same way as for the S 1 x S 2 topology. 



VI. GOWDY MODELS AS SCALAR FIELD THEORIES IN 2+1 CURVED 

BACKGROUND 



The purpose of this section is to reinterpret the reduced models presented in the previ- 
ous sections as certain simple massless scalar field theories in conformally stationary back- 
grounds. We will show how the metrics obtained after the specific gauge fixing and depa- 
rameterization used in the previous sections can be employed to reinterpret the meaning 
(and solution) of the field equations for each topology. This will allow us to use well-known 
techniques of quantum field theory in curved backgrounds to quantize these systems 3(3 . 

Let us start by giving a simple way to solve equations (12.21) 



C a (j)i = 0. 



(6.1) 

(6.2) 
(6.3) 



If a specific solution (g a b,4>i,4>2) is known it is possible to decouple (16. ip and (16. 2116. 3p 
because, when (I6.3P is satisfied, we have the equivalence 



g ab V a V b 



& g ab V a W 



b<Pi 



0. 



The idea is then to solve the last equation in the background g a b and then equation (16. ip 
just gives integrability conditions allowing us to recover g a b. We will discuss next the specific 
form of gab for each of the spatial topologies considered in the paper. 

• Background metric for T 3 

In this case the form of the metric g a b found after the deparameterization is 

p2 t 2 



^a6 = e 7 (-(dt)a(dt)6 + (d0) o (d0), 



b + 



2tt 



-(da) a (da) b 
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defined on (0, oo) x T 2 . A possible (non unique) choice for (<5W>, 0i, 02) is 

9ab = t 2 (-(dt) a (dt) a +(d^) a (d^) 6 + (d ( 7) a (da) ;) ) 

01 = logt 

02 = 

where it is important to notice that even though g ab is not stationary it is conformal 
to a (flat) stationary metric on (0, oo) x T 2 . 

Background metric for S 1 x S 2 

After deparameterization we get now 

g ab = e 7 ( - (d*) a (df ) 6 + (d0) o (d0) 6 ) + ^ sin 2 t sin 2 (da) a (da) 6 
defined on (0, 7r) x § 2 . In this case a convenient choice for (g a b, 4>i, 4>2) is 

g ab = sin 2 t( - (dt) (dt) a + (d9) a (d9) b + sm 2 9(da) a (da) b ^ 

01 = logsin(t/2) - logcos(t/2) 

2 = 0. 

Again this metric is not stationary but it is equal to a time dependent conformal factor 
times the Einstein static metric on (0, it) x § 2 . 

Background metric for S 3 

Finally we have now 

g ab = co S 2 (9/2)e^ - (dt) a (dt) b + (d9) a (d9) b ^ + ^ sin 2 tsin 2 9(da) a (da) b 

defined on (0, tt) x D where D denotes the open disk introduced in the previous section. 
In this possible choice of (g ab , 0i, 2 ) is 

g ab = cos\9/2)e^ - (dt) a (dt) b + (d9) a (d9)^j + sin 2 tsin 2 9(da) a (da) b 
fa = cos6>costlog(tan(t/2)) + cos6> + log(cos 2 (6>/2)) + log(2sint) 

02 = 



where 

sin 2 # 



7 = 



(sin 2 tlog 2 (tant/2) - 2 cost log(tant/2) - l) + log (sin 2 t) 



4 

cost log(tan(t/2)) + cos#costlog(tan(t/2)) + cos6> — I. 



It is important to realize that the concrete functional form of 7 is irrelevant because, 
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whenever £ CT 0j = 0, we have the following equivalence in (0, it) x (§> 2 — {9 = 7r}) 

g ab V a Vb& = & g^VaV^i = 

with 

g ab = sin 2 t( - (dt) (dt) 6 + (d6) a (d6) b + sin 2 9(da) a (da) b ) . 

Notice that the metric g ab is the one that we found for (0, 7r) x § 2 restricted to the 
manifold (0, n) x D obtained by removing a point from the sphere. 

It is important to point out that <pi cannot be extended to the boundary of the disk, 
parameterized as 9 = ir, because (15.111) forces 4>x to behave as log(cos 2 (6'/2)) for 9 — > ir. 
However if we split 0i as 0i = <^ mg + 0^ eg with <^ mg = log(cos 2 (6'/2)) + log(2sint), 
satisfying 

g ab V a V b 0f s = , 

we guarantee that the degrees of freedom contained in 0^ eg still satisfy g ab V ' a V 'b0i Gg = 
(just the same equation as the matter field 02) and can be extended to (0, n) x S 2 . 
Notice that the role of the scalar fields 0^ eg and 02, both well behaved on (0,7r) x S> 2 , 
is symmetric just as in the description of the previous topologies. 



It is important to notice that the scalar field dynamics generated by the time dependent 
Hamiltonians that we have obtained in the previous sections corresponds exactly to the one 
defined by the Klein- Gordon equations on the backgrounds given by g ab . 

To end this section we want to point out that there are certain obstructions to the unitary 
implementation of quantum dynamics for these systems. Specifically, it can be shown that 
it is impossible to find a Fock space representation in which time evolution is unitarily 
implementable [l(J 3). The solution to this problem for the torus case relies on certain 



field redefinitions involving functions of time [7|, [8] . These can be understood in the present 
scheme as coming from the time dependent conformal factors appearing in g ab (or g a b)- As 
we will show in a forthcoming paper, the solution to the unitarity problem for the topologies 
considered here relies on field redefinitions involving precisely the conformal factors shown 
above. Indeed, by performing a redefinition of the scalar fields at the Lagrangian level, such 
that the conformal factor relating both metrics is traded by a time- dependent potential 
term, we expect to find unitary dynamics if this potential is well behaved. In these cases 
the background metric corresponds to a simple, fixed stationary background. 



VII. CONCLUSIONS AND COMMENTS 

We have studied in this paper the Hamiltonian formalism for the compact, linearly polar- 
ized Gowdy models coupled to massless scalar fields. The purpose of the analysis is to have 
a Hamiltonian formulation of the models that can be a starting point for their quantization 
either a la Dirac or by gauge fixing and deparameterization. The results for the T 3 topology 
reproduce the known ones for the gravitational sector and show that the interaction of the 
matter fields occur though the constraints left over by the deparameterization of the system. 
In the other two cases the coupling of matter and gravity degrees of freedom can only be 
seen when the four metric is reconstructed. 
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The description of the S 1 x § 2 and § 3 models requires a careful discussion of the regularity 
conditions that the metric must satisfy in the symmetry axis left over after the Geroch 
reduction performed to describe the systems in 2+1 dimensions. These regularity conditions 
are responsible for the appearance of the so called polar constraints. These can be shown 
to be first class and are necessary conditions to guarantee the differentiability of the other 
constraints present in the models. Of course they must be taken into account in a Dirac 
quantization of the Gowdy models corresponding to these topologies. 

An interesting feature of both the S 1 x § 2 and § 3 cases is the fact that after the deparame- 
terization introduced in the paper (which is a straightforward generalization of the ones used 
in the literature for the familiar T 3 case) there are no constraints left so that the system can 
be completely described by the time dependent Hamiltonians that we have found. This is in 
contrast with the situation for the 3-torus where in addition to the dynamics generated by 
the time dependent Hamiltonian there is an additional constraint in the system that must 
be appropriately taken into account. 

A somewhat surprising fact is the possibility to describe both the S 1 x § 2 and S 3 models by 
using smooth fields on the sphere § 2 . An interesting possibility that may teach us something 
in the case of S 3 is to use a Hopf fibration to perform the Geroch reduction to get a 2+1 
dimensional description. This may be the subject of future work. 

The dynamics of the global modes for the different spatial topologies is easy to obtain 
but there are_ significant differences depending on the topologies. Whereas in the T 3 case 
the value of Q and P are just constants in the other cases P is constant but Q is a function 
of time. 

We have been able to understand in very simple terms the appearance of both initial and 
final singularities in the spacetime metrics that solve the Einstein-Klein Gordon equations 
for these models from the point of view of the phase space description of the dynamics, 
in particular after the deparameterization process that we have followed. As we have seen 
there are natural variables with very simple gauge transformations ( "time dynamics" ) that 
suggest canonical transformations that lead to the time dependent Hamiltonians describing 
the dynamics and explicitly show how the singularities appear. In the S 1 x § 2 and S 3 
topologies the function sin t in a denominator of the time-dependent Hamiltonian shows that 
both final and initial singularities are present whereas the t denominator in the Hamiltonian 
for the 3-torus shows that only an initial (or final) singularity appears in this case. 
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